Genons, twist defects, and projective non-Abelian braiding statistics 
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It has recently been realized that a general class of non-abelian defects can be created in 
conventional topological states by introducing extrinsic defects, such as lattice dislocations or 
superconductor-ferromagnet domain walls in conventional quantum Hall states or topological in- 
sulators. In this paper, we begin by placing these defects within the broader conceptual scheme of 
extrinsic twist defects associated with symmetries of the topological state. We explicitly study sev- 
eral classes of examples, including Z2 and Z3 twist defects, where the topological state with N twist 
defects can be mapped to a topological state without twist defects on a genus g oc TV surface. To 
emphasize this connection we refer to the twist defects as genons. We develop methods to compute 
the projective non-abelian braiding statistics of the genons, and we find the braiding is given by 
adiabatic modular transformations, or Dehn twists, of the topological state on the effective genus 
g surface. We study the relation between this projective braiding statistics and the ordinary non- 
abelian braiding statistics obtained when the genons become deconfined, finite-energy excitations. 
We find that the braiding is generally different, in contrast to the Majorana case, which opens the 
possibility for fundamentally novel behavior. We find situations where the genons have quantum 
dimension 2 and can be used for universal topological quantum computing (TQC), while the host 
topological state is by itself non-universal for TQC. 
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I. INTRODUCTION 

Some of the most important discoveries in condensed 
matter physics over the last few decades have been about 
topological states of matter ii~— Topological states form a 
class of quantum many-body states that are distinguished 
by principles of topology rather than symmetry, which 
implies that they have physical properties that are ro- 
bust under arbitrary local perturbations and that are not 
associated with any local order parameter. Well known 
examples of such topological properties include the exis- 
tence of gapless robust edge states, topology-dependent 
ground state degeneracies, fractional charge and statis- 
tics of topological quasi-particles, and quantized response 
properties. For some topological states, known as non- 
Abelian states^ there are a finite number of degenerate 
lowest energy states for any given configuration of topo- 
logical quasi-particles. The degeneracy between these 
states is topologically protected, being robust against any 
local perturbation. An adiabatic motion of the quasi- 
particles in real space therefore leads to a nontrivial uni- 
tary operation in the space of degenerate ground states, 
which is described by the non-Abelian statistics of the 
quasiparticles. Such a state has been proposed to be 
used, if realized, for intrinsically fault-tolerant, robust 
quantum information storage and processing: the quan- 
tum information is stored in the topologically degener- 
ate states, and the operations are realized by motion of 
quasiparticles. This idea, known as topological quantum 
computation (TQC),— has been a major driving force for 
the study of topological states of matter. 

The discussion of topological quasiparticles in topolog- 
ical states can be extended to extrinsic defects, which are 
point-like objects that are not intrinsic, finite-energy ex- 
citations of the system, but instead have a long-ranged 
confining interaction with each other£~— Such defects 
can be created and controlled by an external field, and 



they may also carry non-abelian statistics, similar to in- 
trinsic topological quasiparticles. However the overall 
phase of the unitary operation generated by braiding the 
defects is generically non-universal and path dependent, 
due to their long ranged interaction. Therefore the statis- 
tics of the extrinsic defects is only well-defined up to a 
phase, implying that their braiding forms a projective 
representation of the braid group. This possibility is re- 
ferred to as projective non-Abelian statistics^ A simple 
example of a non-trivial extrinsic defect is the vortex in 
a two-dimensional p x +ip y topological superfluid 5 , which 
has Majorana zero modes and non-Abelian statistics sim- 
ilar to the topological quasi-particles in the Moore-Read 
Pfaffian fractional quantum Hall (FQH) stated The 
overall phase of the statistics is undetermined due to the 
logarithmic interaction of the vortices. Recently, similar 
extrinsic defects with Majorana zero modes have been 
discussed in several other physical systems^. 

Since projective non-abelian statistics are not de- 
scribed by the same mathematical framework 15 as in- 
trinsic topological quasiparticles, there may potentially 
be a world of possibilities that would be inconsistent 
for ordinary non-abelian statistics. An important direc- 
tion therefore is to develop the theory of projective non- 
abelian braiding statistics and to investigate the novel 
possibilities. 

A second important direction is to explore ways of 
obtaining exotic non-Abelian defects, beyond Majorana 
fermions, in simple, experimentally achievable settings. 
This is motivated in part by the well-known deficiency 
of Majorana fermions to yield a universal gate set for 
TQC. 1 — Recently, extrinsic defects with projective non- 
Abelian statistics beyond Majorana fermions have been 
proposed in several systems. In fractional Chern insula- 
tors (i.e. FQH states in lattice models without an exter- 
nal magnetic field)^£r— , a class of states called topologi- 
cal nematic states can be realized when the Chern num- 
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ber of the partially occupied energy band is larger than 
l£ Topological nematic states are topologically equiva- 
lent to conventional multi-layer FQH states on regular 
lattices, but lattice dislocations effectively change the 
topology and introduce "worm holes" between the two 
layers, which leads to the non-Abelian statistics of the 
dislocations. Similar non-Abelian lattice dislocations can 
also be realized in rotor modelsi£i2ii£. A different 
type of extrinsic defect has been proposed on the bound- 
ary of FQH or fractional quantum spin Hall (FQSH) 
states, at domain walls between regions in which the edge 
states obtain different mass termsi22r— 

In this paper, we develop a more systematic under- 
standing of the projective statistics of these extrinsic de- 
fects and the relation between different types of extrinsic 
defects. We study a wide class of defects that obey pro- 
jective non-abelian statistics, and that can be interpreted 
as twist defects associated with symmetries of the topo- 
logical state. The examples discussed in the last para- 
graph can all be understood as twist defects. For ex- 
ample, a Chern number 2 topological nematic state is 
mapped to a bilayer FQH state, with a Z2 symmetry 
associated with exchanging the two layers. The lattice 
dislocation is a Z2 twist defect in the sense that a quasi- 
particle going around the dislocation once will be acted 
upon by the Zi operation of exchanging the two lay- 
ers. In the simplest case, a quasi-particle in one layer 
will end up in the other layer upon winding around the 
dislocation. 9,26 Such a Z2 twist defect can also be gener- 
alized to any other Abelian or non-Abelian bilayer topo- 
logical states. Similarly, a Zm topological state has a Z2 
symmetry associated with electric-magnetic duality, al- 
lowing for the possibility of twist defects that exchange 
electric and magnetic quasiparticles as they encircle the 
defect ^ 10 ' 26 Other examples that we will study in this 
paper include the twist defects associated with a certain 
particle-hole symmetry in FQH states, which are equiv- 
alent to the edge defects in FQH/FQSH systems studied 
recentl y 10 ' 22 ~ 25 i 27 , and the Z% twist defects that can ap- 
pear in triple-layer FQH or topological nematic states. 

We develop several complementary methods of com- 
puting the projective non-abelian braiding statistics of 
these twist defects. We find that in the class of examples 
that we study, the state with N twist defects can always 
be mapped - in a certain precise sense - to a topolog- 
ical state without twist defects, but on a genus g oc N 
surface. Then, we find that the braiding of the twist de- 
fects realizes adiabatic modular transformations, or Dehn 
twists, of the topological state on the high genus surface. 
This provides a physical way to implement elements of 
the mapping class group of a topological state on a high 
genus surface. Since adding twist defects effectively in- 
creases the genus, we refer to them as genons^ 

One particularly interesting example is given by genons 
in an Isingxlsing topological state. The Ising topologi- 
cal state is a simple non-Abelian topological state with 
three topological quasiparticles, which can be realized in 
Kitaev's honeycomb lattice model 2 ^.. By Isingxlsing, we 



are referring to a bilayer state with each layer correspond- 
ing to an Ising theory. It is known that braiding particles 
in the Ising theory is not sufficient for universal topologi- 
cal quantum computation. 15 However, we will show that 
utilizing the braiding of twist defects in an Isingxlsing 
theory can make the state universal, since it realizes Dchn 
twists of a single Ising theory on high genus surface a 29 ' 30 . 
We will also demonstrate that the average degree of free- 
dom, i.e., the "quantum dimension", of each twist defect 
in this theory is d — 2. This provides an interesting ex- 
ample where a defect with integer quantum dimension 
can allow for universal TQC while the host topological 
state is by itself non-universal. 

The genons are confined in the sense that there is a 
long-ranged confining potential that grows with their sep- 
aration. One way to deconfine them is to gauge the sym- 
metry they are associated with.— This leads to a new 
class of topological states ) 26 ' 31 i 32 where the defects are 
now intrinsic non-abelian quasiparticles. In this paper, 
we discuss the relation of the braid matrices between the 
cases where the genons are confined, and cases where they 
are deconfined by gauging the associated symmetry. In 
the case where the genons have qauntum dimension V% 
they can be interpreted as Majorana fermions, and gaug- 
ing the symmetry does not change the braid matrix, but 
only makes the overall phase well-defined and universal. 
In the more general situations, we find that gauging the 
symmetry changes the dimension of the braid matrix, so 
that strictly speaking the braiding is different, though 
closely related. In the single-component case, we find 
that gauging the symmetry can even change the quan- 
tum dimension of the non-abelian defects. These provide 
simple examples where projective non-abelian braiding 
statistics can give braiding that is inequivalent to ordi- 
nary non-abelian statistics. 

The rest of this paper is organized as follows. In Sec. 
ITTl we introduce the notion of a twist defect in a topolog- 
ical state and we discuss several examples, possible phys- 
ical realizations, and we give a brief discussion of projec- 
tive non-abelian statistics. In Sec. IIII1 we study in de- 
tail the braiding of Z2 twist defects in "two-component" 
states, which can described by U(l) x U(l) CS theory and 
which include double layer FQH states and Zjy topologi- 
cal states. We discuss the sense in which the twist defects 
can be thought of as "genons" (subsequently, we will use 
"twist defect" and "genon" interchangeably in this pa- 
per). In Sec. IIV[ we study Z2 genons in single-component 
states, which are described by U(1)n CS theory, and we 
discuss the close relation between these genons and the 
ones in the two-component case. In Sec. |V1 we study Z2 
genons in decoupled, double-layer non-abelian states, and 
discuss the possibility of universal TQC with genons. In 
Sec. IVI1 we study Z3 genons in three-component FQH 
states, which are described by U(l) x U(l) x U(l) CS 
theory; this provides the first example of projective non- 
abelian braiding beyond the Z2 genon case. In Section 
IVIII we discuss topological states that are obtained when 
the symmetry associated with the genons are gauged, and 
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we discuss the relation between the braiding of the genons 
before and after the symmetry is gauged. We conclude 
with a discussion in Section IVIIII 

II. TWIST DEFECTS IN TOPOLOGICAL 
STATES 

A. General definition of twist defect 

A topologically ordered phased is generally charac- 
terized by a set of topologically non-trivial quasiparti- 
cles, {ji}, for i = 1, • • • , N qv , where N qp is the number 
of quasiparticles. Below, we will briefly sketch the topo- 
logical properties of the quasiparticles. First, when two 
quasiparticles are observed from far away, they behave 
like a superposition of single quasiparticle states. This is 
described by the fusion rules 7; x 7^ = J2 k N^k- Sec- 
ondly, when two quasiparticles ji , jj wind around each 

other, a phase e ie ^ is obtained, which depends on the 
fusion channel k. 9% is referred to as the braid statistics 
of the quasiparticles. When a particle is spinned around 
itself by 2tt, it generically gains a nontrivial phase e l6i . 
6i = for bosons and tt for fermions, and in general it 
can take any value between [0, 2n). The braiding, fusion 
rules and spins need to satisfy some consistency condi- 
tions but we will not review them here.— Mathematically, 
a topologically ordered state is characterized by a unitary 
modular tensor category (TJMTC) J£>22££ 

It is possible for a topological phase to have a discrete 
symmetry g which maps a quasiparticle ji to another 
particle 7 sW , 

7i -> lg(i), £=!,••• ,N qp (1) 

while preserving all topological properties such as fusion 
rules, braiding and spins. All such symmetries form the 
group of automorphisms of the UMTC. Table U summa- 
rizes some examples. For instance, a Zn topological state 
has N 2 quasiparticles, which can be labelled as (a, b), for 
a, b = 0, • • • , N — 1. The (a, 0) particles are the electric 
particles, while the (0, a) particles are the magnetic ones. 
This state has a Z2 x Z2 symmetry. One of the Zi sym- 
metries is the electric- magnetic duality (a,b) — > {b,a), 
ie exchanging electric and magnetic particles. The other 
Z2 is associated with taking (a, b) —> (N — a, N — b), 
which takes the electric and magnetic particles to their 
conjugates. 

If charge conservation is broken, a 1/ fc-Laughlin FQH 
state also has a Z2 symmetry, associated with exchanging 
quasiparticles and quasiholes. This is beacuse the quasi- 
particles and quasiholes have the same fractional statis- 
tics and yet are topologically distinct quasiparticles. 

In a bilayer FQH state, there can be a Zi symmetry 
associated with exchanging layers. In an A^-layer FQH 
state the symmetry of permuting the layers is Sn, which 
may be broken to a smaller subgroup, such as Zjy. 

Given a topological phase with such a symmetry, a 
twist defect is an extrinsic defect of the system, labelled 




FIG. 1: Worldline of quasiparticle 7, and twist defect, la- 
belled by g. Braiding ji around the twist defect changes the 
quasiparticle to J g (i)- 



by the symmetry g, such that a quasiparticle 7* gets 
transformed to 7 s (j) as it braids around the twist defect 
(see Fig. [1]). Such twist defects generally have a non- 
trivial quantum dimension and lead to topological degen- 
eracies. Therefore, every automorphism of the UMTC 
is associated with a twist defect that realizes projective 
non-abelian statistics. 

The twist defects discussed above are related to recent 
discussions in the mathematical physics literature^i^r— 
where the mathematical theory of twist defects and 
boundaries between topological states is currently under 
investigation. 

There are also twist defects that cannot be included in 
the definition above — those that act trivially on the topo- 
logically distinct quasiparticles, but have a non-trivial 
action on local operators. In these situations, it is possi- 
ble but not guaranteed that the twist defect will have a 
non-trivial quantum dimension. One example is a system 
with charge conservation, where the twist defect has the 
effect of sending a quasiparticle to its conjugate. In an 
IQH state, such defects can trap Majorana zero modes, 
even though the IQH state does not have any topolog- 
ically non-trivial quasiparticles. A second example is a 
superconductor, with the Z2 symmetry that takes the 
electron c — > — c. The vortex in p + ip topological super- 
conductors, which traps a Majorana zero mode, can then 
be interpreted as a twist defect that acts on the fermions 
by this Z2 symmetry as they encircle the vortex. 



B. Physical Realizations 

Recently, twist defects have been shown to occur in a 
number of different physical realizations. These include 
dislocations in topological nematic states realized in FCI 
with higher Chcrn number C > l 9 -, at certain junc- 
tions in gated bilayer FQH states^, dislocations in ex- 
actly solvable Zm topologically ordered models£&i£, and 
superconductor-ferromagnet domain walls at the edge of 
2+ld TIs and FTIs l 2fc242I 

The topological nematic state realization of the twist 
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Examples of Symmetries of Topological States 

Topological states Symmetries Transformation of quasiparticles 

„ Electric-magnetic duality Zi (a, b) — > (b, a) 

Zj jv states 

Particle-hole symmetry Zi (a, b) — > (jV — a, N — b) 

iV-layer FQH states Layer permutation Sn (oi, 02, ••, cln) (iPi, ap 2 , ap N ) 

1/fc-Laughlin FQH state Particle-hole symmetry Z2 a — > (k — a). 



TABLE I: Some examples of topological phases and their symmetries. Zn topological states have Z2 x Zi symmetry associated 
with both the electric-magnetic duality ( (a, b) — > (b, a) ) and the particle-hole transformation taking quasiparticles to their 
conjugates: (a,b) — > (iV — a,N — b), where a,b — 0, ■ ■ ■ , N — 1. iV-layer FQH states can have a symmetry associating with 
permuting layers. 1/fc-Laughlin states can have a particle- hole symmetry associated with taking a quasiparticle to its conjugate, 
if charge conservation is broken. 



defects^ is based on the Wannier state representation of 
FCl2£. This representation established a mapping be- 
tween one-dimensional Wannier states of an FCI system 
and the Landau level wavefunctions in the Landau gauge 
in an ordinary FQH system. By such a mapping, a Chern 
number C = 1 band is mapped to a Landau level. The 
validity of the Wannier state representation approach has 
been confirmed in C = 1 states with filling 1/2 and 
l/322i&. FCI with C > 1 bands have also been studied 
in analytic and numerical work o 9 i 41 ~ — . In the Wannier 
state representation, a band with Chern number C > 1 
is mapped to a C-layer quantum Hall system^. In con- 
trast to ordinary multi-layer quantum Hall systems, the 
different effective layers in this system are related by lat- 
tice translations, which implies that lattice dislocations 
can act as twist defects: as a reference point encircles the 
dislocation, it is translated by the Burgers vector of the 
dislocation and may effectively belong to a different layer 
of the effective FQH system. 

The realization of twist defects as lattice dislocations 
in the Zn models is similar. 7,10 In some lattice models of 
the Zn states, the electric particles {Zn charges) belong 
to the even sublattice, while the magnetic particles (Zn 
vortices) belong to the odd sublattice. A lattice disloca- 
tion can create a situation where an excitation that encir- 
cles it starts in the even sublattice and ends in the odd 
sublattice, which therefore implies that the dislocation 
can act as a twist defect whereby electric and magnetic 
particles are exchanged upon encircling it. 

As we will explain in subsequent sections, the 
superconductor-ferromagnet domain walls at junctions of 
IQH/FQH states can be viewed as twist defects in several 
different ways: they can be viewed as twist defects asso- 
ciated with a Z2 particle-hole symmetry of the Abelian 
QH states, or they can be mapped onto double layer FQH 
states, in which case they are associated with the Z2 layer 
exchange symmetry. 



C. Projective non-abelian braiding statistics 

The twist defects are extrinsic defects of the topologi- 
cal phase, not finite-energy, intrinsic excitations. There- 



fore, separating twist defects will typically cost an energy 
that grows with the distance between them, either log- 
arithmically or linearly, depending on whether the asso- 
ciated symmetry is a continuous symmetry of the state. 
This is similar to the vortices in topological superfluids 
(and thin-film topological superconductors), which also 
have a logarithmic energy cost associated with separat- 
ing vortex/anti- vortex pairs. In these cases, we say that 
the defects are confined . 

Since the defects can have a non-trivial quantum di- 
mension, braiding them can lead to non-abelian braid- 
ing statistics. However, since well-separated twist de- 
fects still have an energy cost that grows with the sepa- 
ration, the overall phase of the braiding statistics is not 
well-defined. To see this, observe that simply moving 
one quasiparticle in a small circle without encircling any 
other quasiparticle still accumulates a phase from the 
non-negligible change in energy of the state during the 
process. Since the overall phase of the statistics is not 
well-defined, these twist defects form a projective repre- 
sentation of the braid group. In contrast, in a conven- 
tional non-abelian state, the non-abelian quasiparticles 
are finite energy, deconfined excitations, and the overall 
phase of the braiding is topological. Therefore, the braid- 
ing of non-abelian quasiparticles in a true non-abelian 
state forms a linear representation of the braid group. 

The quasiparticles of a true non-abelian state are sub- 
ject to the mathematical constraints of a UMTC. How- 
ever, the twist defects, since they form a projective rep- 
resentation of the braid group and are not intrinsic quasi- 
particles, are described by a different mathematical the- 
ory (see eg Ref. [l2l ) , and in principle are not constrained 
in the same way. This opens the possibility of fundamen- 
tally novel behavior. 



III. Z 2 TWIST DEFECTS IN 
TWO-COMPONENT ABELIAN STATES 

In this section, we concentrate on the properties of Z2 
twist defects in Abelian states that can be described by 
U(l) x U(l) CS theory. This includes two-component 
FQH states^ and Zn topological states. The Lagrangian 



FIG. 2: A pair of twist defects induces two distinct non- 
contractible loops, labelled a and b. For twist defects in 
(mmO) states, these loops effectively cross only once, leading 
to a magnetic algebra for the quasiparticle loop operators. 

is given in terms of two U(l) gauge fields, a and a: 

7TL I 

C = — (ada + ada) + — (ada + dda), (2) 

47T 47T 

where m and I are any integers and ada = a^d^ax. 
To describe Zn topological states, we set m = and 
I = N above. This theory encodes the topological prop- 
erties of these topological states. As we review below, the 
Z2 twist defects in these states a re non-A belian defects, 
carrying a quantum dimension \J\m — ^ 

A. Bulk geometrical picture 

The twist defects introduced above are point defects: 
far away, no local operator can distinguish their pres- 
ence. However, for the purpose of understanding the be- 
havior of the twist defects, it is helpful to imagine that 
the twist occurs along a single branch cut that connects 
them. This is similar to supposing that the phase wind- 
ing of a vortex in a superconductor is all localized to a 
single cut connecting a vortex/anti- vortex pair. 

The twist defects introduce new non-contractible loops 
into the system, along which the quasiparticles can prop- 
agate. For example, with two pairs of Z2 twist defects 
on a sphere, there are two distinct non-contractiblc loops 
(Fig. [2]). The key feature of the twist defects is that their 
presence yields a non-trivial algebra for quasiparticle loop 
operators corresponding to these non-contractible loops. 
For example, consider the case of two decoupled 1/m 
Laughlin FQH layers (denoted (mmO) states), where the 
twist defect exchanges the layers. The two loops a and b 
effectively cross only once, because a quasiparticle goes 
from one layer to another as it passes through the branch 
cut, and the quasiparticles in two different layers have 
trivial mutual braid statistics. This leads to the mag- 
netic algebra 

W(a)W(b) = W(b)W(a)e 27rl/m , (3) 

where W(C) is the operator that tunnels a Laughlin 
quasiparticle around the loop C . The explicit expression 
of W(C) is given later in Eq. ([7]). The ground states 
form an irreducible representation of this quasiparticle 
loop algebra, which in this example is |rn|-dimensional. 



FIG. 3: n pairs of twist defects induces 2(n — I) distinct 
non-contractible loops, at, bi, for i = 1, • • • , n — 1. For twist 
defects in the (mmO) states, the quasiparticle loop operators 
give rise to n — 1 copies of the magnetic algebra, leading to 
|m| n_1 topologically degenerate states. 




FIG. 4: The (7(f) x (7(1) CS theory with n pairs of disloca- 
tions on a sphere can be mapped to a (7(1) CS theory on a 
genus g = n — 1 surface, M„-i. M n -i consists of two copies 
of the original space. A new (7(1) gauge field, c, is defined on 
Mn-i, such that c = a on the top half, and c = —a on the 
bottom half of M„_i.- 



Therefore we can conclude that the twist defects have a 
non-trivial quantum dimension. 

Proceeding with the example of the (mmO) states, 
with n pairs of twist defects on a sphere, we can de- 
fine 2(n — 1) non-contractible loops (see Fig. [3]), a*, b, 
for i = 1, • • • , n — 1, such that and bj effectively cross 
exactly once if i = j, and do not cross otherwise. 

Using these non-contractible loops, the quasiparticle 
loop algebra now becomes 

W(cn)W(bj) = W{b j )W{a i )e 5 ^' Ki / m . (4) 

Thus we have n — 1 copies of the magnetic algebra, and 
therefore a finite dimensional irreducible representation 
of dimension |m|™ -1 . This shows that the quantum di- 
mension of each twist defect in the (mmO) states is -\/|m|. 

It is possible to derive the above results more con- 
cretely by starting with the field theory @. In the 
(7(1) x (7(1) CS theory, a pair of Z 2 twist defects can 
be modelled as a pair of points, connected by a branch 
cut 7 such that at the branch cut, the two gauge fields 
a and a obey twisted boundary conditions. Defining 
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A = { a ^ | , this means 
\0 a ' 



lirn A(p) 



lim a x A(p)a x , 



(5) 



for every point po on 7^ The limit p — > p^ ' means that 
the limit is taken approaching one particular side (or the 
other) of 7. This ensures that quasiparticles encircling 
a twist defect get transformed by the Z2 action, o~ x — 

(: :) 

In the presence of n > 1 pairs of twist defects, it is pos- 
sible to consider a single U(l) gauge field c on a doubled 
space, M n _i, where M n _i is a genus g = n — 1 surface 
(see Fig. @|. As shown in Fig. [U c = a for points in the 
top half of M n -i and c = —a for points in the bottom 
half of M„_i. The original £7(1) x U(l) CS theory © 
can then be shown to equivalent to a U(l) m -i CS theory 



is then given by non-abelian adiabatic Berry phases as- 
sociated with Dehn twists, which were computed in Rcf. 
0for U(l) CS theory. Note that as discussed in Ref. 0, 
these non-abelian adiabatic Berry phases are well-defined 
only up to the overall phase of the matrix, reflecting the 
fact that the non-abelian braiding statistics are projec- 
tive. 

Using the relation to Dehn twists, we can now imme- 
diately obtain the braid matrices. Consider two pairs of 
twist defects on a sphere, so that we are mapped to a 
U(l) m -[ CS theory on a torus. Let us consider a basis 
that consists of wrapping the quasiparticles around the b 
direction of the torus: 



|n), n = 0, • • • , \m — l\ — 1. 



(10) 



This corresponds to diagonalizing W(a), while W(b) acts 
as a raising operator: 



W(b)\n) = \n+l mod \m-l\), 
W(a)\n) = e 2m/{m - l) \n). 



(11) 



47T 



cdc. 



(6) 



M„_i 



The quasiparticle loop operators W(C) defined earlier in 
the context of the (mmO) states are written in the field 
theory as 



W(C) =Ve i fc c -< a , 



(7) 



Such a theory has a ground state degeneracy of \m— Z| n_1 , 
which shows that th e Z2 twist defects have a quantum 
dimension yfro — 

If the n pairs of twist defects were placed on a genus g 
surface instead of a sphere, then the ground state degen- 
eracy would be given by |m 2 — l 2 \ 9 x \m — l\ n ~ l , where 
the factor |m 2 — l 2 \ 9 comes from the "bare" degeneracy 
of the manifold on which the twist defects are placed. 

The explicit mapping to the high genus surface is par- 
ticularly useful for understanding the non-abelian braid- 
ing statistics of the Z2 twist defects. First, we exchange 
the twist defects on the sphere, and then we can visual- 
ize how the non-contractible loops transform under the 
exchange. From Fig. [S^, andJSJj, we can see that under a 
counterclockwise exchange B12 of the defects labelled 1 
and 2, the loops transform as 



B\2 



(8) 



From Fig. [SJ; and [SJi, we see that under clockwise ex- 
change of the defects labelled 2 and 3, 



D. 



(9) 



Therefore, we see explicitly that the braiding of the twist 
defects correspond to adiabatic modular transformations 
(Dehn twists) in the effective genus g surface. The adia- 
batic Berry phase associated with braiding twist defects 



Dehn twists along a non-contractible loop C will be de- 
noted by Uc- Based on the previous discussion, 



B12 — U a , B 2a — Ub- 



(12) 



It was found^ that U a are diagonal in the above basis, 
and given by 



UJn) 



^i9giirn / (m-l) + iirn(m-l) i \ 



(13) 



where the overall phase e depends on details of the path, 
not only on its topology. In order to compute we 
observe that 



Ut = S^U a S, 



(14) 



where S is the modular S matrix, which exchanges the a 
and b cycles: 



S : a —¥ 6, b —¥ —a. 



(15) 



The adiabatic modular transformation associated with S 
is given by^ 



S, 



1 



a/3 



a -2-7ria/S/(-m-i) 



(16) 



Therefore: 



1 1 /3=0 



iTr(/3 2 +2P(a-S))/(m-l) 



(17) 



where the overall phase e 1 ^ again depends on non- 
universal details of the path. 

Therefore, we find that in the presence of the twist 
defects, there is a group of projective non-abelian braid- 
ing statistics associated with the braiding of the defects, 
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FIG. 5: (a) A loop that encloses the twist defects labelled 2 and 3 is mapped to the b cycle of the torus, (b) Effect of a 
counterclockwise exchange of 1 and 2. By following the effect on the non-contractible loop, we see that in terms of the genus g 
surface, it has the effect of a Dehn twist along an a cycle. Thus, the original b loop becomes the loop b + a after counterclockwise 
exchange of 1 and 2. (c) A loop that encloses a pair of twist defects connected by a branch cut is an a-cycle of the genus g 
surface, (d) Effect of a clockwise exchange of 2 and 3. We see that the a loop gets mapped to the a + b loop. Thus the clockwise 
exchange of 2 and 3 is equivalent to a Dehn twist along the b loop. 



along with the braiding of the quasiparticles around the 
defects. These are generated by the unitary operators 

{B iti+l ,W{ai),W{bi)}. (18) 

The generalization of the braid matrices to n > 2 pairs 
of twist defects is straightforward, since for the U(l) m -i 
CS theory, the states associated with each handle of the 
genus g — n — 1 surface are independent. 

We see that the effect of the n twist defects is to intro- 
duce a set of non-trivial quasiparticle loop operators that 
satisfy the same algebra as the quasiparticle loop oper- 
ators of a U(\) m ~i CS theory on a genus n — 1 surface, 
and that braiding of the twist defects is directly related 
to elements of the mapping class group of this theory 
on the high genus surface. In order to emphasize the 



connection between the twist defects and the high genus 
surface, we refer to the twist defects as genons. All of the 
twist defects considered in this paper have this property, 
and therefore in this paper we will use "twist defects" 
and "genons" interchangeably. 



B. An example: Majorana braiding as Dehn twist 
of 1/2 Laughlin state 

The simplest case where genons are non-abelian is for 
the case \m — l\ — 2, in which case they have quantum 
dimension We therefore expect that a Majorana zero 
mode is localized at the defect, and the braiding statistics 
should be associated with the braiding of Ising anyons. 



Let us consider the case where we have n = 2 pairs 
of genons on a sphere, which maps us to a genus g = 1 
surface. In this case, there are \m — l\ 9 = 2 ground states. 
The braiding of the genons, given by the Dehn twists are: 



U a = e 





(19) 



Observe that U a is, up to an overall phase, the same as 
the braid matrix of Ising anyons!^ In other words, Dehn 
twists in the U{V)2 CS theory are equivalent to projective 
Ising braiding statistics, due to this profound relation to 
Z2 genons. 

In addition to the braiding of the genons, we can also 
braid quasiparticles around genons to get: 



W(a) 







W(b) 




(20) 



C. 1 + 1D edge CFT picture 

It is helpful to derive the above results in a differ- 
ent way using the one-dimensional chiral Luttinger liquid 
edge theory. This will give us a different perspective on 
how to compute the topological degeneracy and braiding 
statistics, it will connect to the more standard zero mode 
analysis used in other contexts^ and it will help give a 
different protocol for carrying out the braiding. The lat- 
ter will be useful depending on the physical realization 
of the genons. First, we align all of the defects along a 
single line, and we cut the system along the line. Then, 
we have two counter-propagating chiral Luttinger liquid 
theories, with the left-moving one localized on one edge, 
and the right-moving one on the other edge. The edge 
theory for U(l) x U(l) CS theory described by a generic 
if-matrix is given by the actior 



,45,46 



s, 



edge 



1 

An 



dxdt[Kud t (f> L id x <p L j - Vud x <t> L id x (j>Lj\, 

(21) 



where in our case K 




t>Li denotes left- 



moving chiral bosons for I = 1, • • • , dim KM- Here and 
below the repeated indices J, J are summed. The field 
</>lj is a compact boson field with radius R = 1: 

4>LX ~ <t>Li + 2tt. (22) 

Quantizing the theory in momentum space yields* 

[d K 4>Li(x),4>Lj(y)] = -i2nKjjS(x - y). (23) 

Integrating the above equation gives: 



■'LI 



(x), <f> LJ {y)} = -i-KKj}sgn(x - y) (24) 




(b) 



H-jr u T J 




1 












1 fr" 




* t ! ; r . 






-4 PU 











FIG. 6: The edge state understanding of the topological de- 
generacy, (a) The genons are oriented along a single line, 
and then the system is cut along the line, yielding gapless 
counterpropagating edge states along the line. The origi- 
nal topological state is obtained from gluing the the system 
back together by turning on appropriate inter-edge tunneling 
terms, (b) Depiction of the two branches (red and blue) of 
counter-propagating edge excitations. The arrows between 
the edge states indicate the kinds of electron tunneling terms 
that are added. Away from the genons, in the A regions, the 
usual electron tunneling terms involving tunneling between 
the same layers, ' eRI i S/eLi + H.c, are added. In the regions 
including the branch cuts separating the genons, twisted tun- 
neling terms are added: ^ eR1 9 e L2 + ^\r2^<^li + H.c. 



The electric charge density associated with <f>Lj is given 
by 



PLI = —d x <j) L I, 
ITT 



(25) 



and the 7th electron operator is described by the vertex 
operator 



eLI 



JKi 



(26) 



Note that normal ordering will be left implicit (ie 
giKu^Lj =■ e iKu<t> L j .y jf we cons ider the FQH state 

on a cylinder, we will have a left-moving chiral theory on 
one edge, and a right-moving chiral theory on the other 
edge. For the right-moving theory, the edge action is 

Sedge = / dxdtl-Kljdt^Rldx^RJ - Vljd x (f>Rld x 4>Rj], 

(27) 

so that 

[<j)Ri{x), (j)Rj{y)] = mKjjsgn{x - y), (28) 
the charge is 



Pri = -z-dx<t>Ri, 

Z7T 



(29) 
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and the electron operator is 

^eRI 



— e —iKij<j>nj 



(30) 



When m is odd, the electron operators are fermionic, 
so we need to ensure that ^ e Ri{x)^ e Lj{y) — 
e mK ' J ^eLj{y)^eRi{x)- This can be done by introduc- 
ing commutation relations: 

[cj) LI (x),<p R j{y)} = ittK~J. 



(31) 



Now suppose we have n pairs of genons, such that going 
around each genon exchanges the two layers. Each pair 
of genons is separated by a branch cut. Let us align all 
the genons, and denote the regions without a branch cut 
as Ai, and the regions with a branch cut as Bi (Fig. [BJV). 
Now imagine cutting the system along this line, introduc- 
ing counterpropagating chiral edge states. The gapped 
system with the genons can be understood by introduc- 
ing different electron tunneling terms in the A and B 
regions (Fig. [BJ3) : 



KL2^eR2 



eRl 

Introducing the variables 



>L1 
t>L2 



»R2 
i>Rl 



H.c if x e Ai 
H.c if x € B t 

(32) 



(33) 



we rewrite (1321) as 



5U t 



J2iCos{Ku<j>j) if x e Ai 
X^cos^/ j0j) if x G Bi 



(34) 



It is helpful to rewrite the above as 
'to + I 



<% = g cos 




cos(^0_) if xeA t 
cos(2^i0_) if xeB t 



(35) 



where <j>± — (f>\ ± fa, and cj>± — (f>\ ± <f>2- One way to 
understand the topological degeneracy was explained in 
Ref. [^. In the absence of the twist defects, there are 
\(m + l)(m — l)\ states, associated with the distinct eigen- 
values of e <Oi±<fc»), which are given by e 2 ™P±/( m± '), for 
p± integers. Physically, e l(t " corresponds to a quasiparti- 
cle tunneling process, where a quasiparticle from the 7th 
layer is annihilated at one edge, tunnels around the torus, 
and is created at the other edge. In the presence of the 
n pairs of twist defects, the eigenvalue of e l( ^+ is globally 
pinned everywhere (see eq. (|3"5|) ) while e 1 ^- can take 
\m — l\ different values in each of the Ai regions. The op- 
erator e 4 ^ 1- ^ 2 ) is unphysical and not gauge-invariant,— 
so we cannot label the states by its eigenvalues in the B 
regions. This yields a total of \(m + l)(m — l) n \ states and 
agrees with the bulk calculation of Sec. MI AI with the 
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FIG. 7: Physical meaning of the zero mode operators at and 
Pi. In subsequent figures, we will use the bottom figure, which 
is more precisely described by the top figure. 



extra factor of |to 2 — l 2 \ due to the fact that in this case 
the defects were placed on a torus instead of a sphere. 

In order to calculate the braiding and to understand 
the connection to zero modes localized at the defects, 
we consider quasiparticle tunneling operators near each 
genon: 



i<t>i(x; 



-i<pi(x Bi ) 



Oi2i-\ = e ~<'e 

_ e i<t>l(x Bi ) e -i<p2(xA i+1 ) 

J4>1 (x Bi ) e ~ i( t>l (x A i+1 ) 



a-2i 



(36) 



for i = 1, • • • , n, where XAi and XB t are the midpoints of 
the Ai and Bi regions, respectively. Physically, these op- 
erators are quasiparticle tunneling operators, projected 
onto the ground state subspace where the fields are con- 
stant within each region. Fig. [7] displays the quasiparti- 
cle tunneling process described by these operators. The 
operators ccj and (3i are zero modes: 



(37) 



where H is the edge Hamiltonian including the tunneling 
terms (f3"2")) . In the ground state subspace, as discussed 
above, the value of e l ^ 1+ ^ is pinned everywhere. Thus: 



OL2i+lft2i+l — lj 



OL2i$2 



(38) 



Therefore, there is a single independent zero-mode oper- 
ator, OLi, which satisfies an algebra: 



z2tt / (m — l) 



(39) 



where k > 0. For |m — l\ = 2, ai are Majorana fermions. 
For | to — l\ > 2, this can be viewed as a generalization of 
the Majorana fermion algebra, which is usually referred 
to as a Zi m _n parafermion algebra4£ 

Alone, the zero mode operators ai are not physical 
gauge invariant operators, as they correspond to the 
quasi-particle motion along an open path. However, 
proper combinations of them are physical gauge invariant 
operators. For example, 

W{(a i+ i) = a 2l +i^ 2i +2, W 2 '(a i+ i) = /3 2l +ia 2 i + 2 (40) 
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FIG. 8: Effect of a clockwise braid on the zero modes a 1 and 
Qf2- (a) Qi — > oi2. (b) Q2 gets transformed to a combina- 
tion of two loop, /3i/3|, an d OL2- Therefore, Q2 — > ct20e\a2, 
where recall j3i — a\. As explained in the text, there is an 
ordering ambiguity between 02 and a\a2 that is fixed with 
an additional constraint. 

describe quasiparticles from layer 1 and 2, respectively, 
tunneling around the pair of genons 2i + 1, 2i + 2. Simi- 
larly, 

Wl(b i )=a 2i p 2i+ i, W%(b i+1 ) = p 2i a 2i+ u (41) 

where a t and 6j are as shown in Fig. [3J Using the algebra 
of these quasiparticle tunneling operators in the edge the- 
ory, we can obtain the algebra of the quasiparticle loop 
operators 

WMWiih) = w^hWM)^ 1 '^- (42) 

Note that from ([35]), it follows that Wi(C) oc W%{C). It 
is useful to note that 

n m—l 1 m—l { 1 \m — I — 1 f ao\ 

a 2i+i = 1> a 2i = > (43) 

so that 

W™- l (a l ) = W™-\b l ) = l. (44) 

The irreducible representation of this algebra contains 
|m — l\ n ~ x states. This is simply the 1+1D edge CFT un- 
derstanding of the bulk geometric construction described 
in Section InTAl' 

Using the zero mode operators cvi, we can develop a 
1+1D CFT understanding of the braiding of the twist 
defects. First, consider how the quasiparticle loop oper- 
ators are transformed under a clockwise exchange of the 
defects 1 and 2, which we denote by the unitary operator 
B\ 2 . The result can be inferred simply by the diagrams; 
we draw the physical processes associated with a,\ and 
«2j and then we consider the effect of the clockwise braid, 
keeping the open ends fixed. From Fig. [5J we can see that 



B\ 2 a\Bi2 = e %v OL 2 and B\ 2 a 2 B\ 2 = e z8 a2a\a2- We can 
fix the relative phase e lv = e lB by using the fact that 
the braiding of 1 and 2 should keep invariant the eigen- 
values of the loop operator W\{a) = a\a2 that encircles 
the pair. We can partially fix the remaining phase by 
observing that the Z2 layer exchange symmetry implies 
BlzfcBu = e*fo and B\ 2 p 2 B l2 = e^hftfa Using 
([35]). we can then determine that e l2ip = e -2W(m-0. 
Thus, we find: 

B\ 2 a!B 12 = e «^-W('»-0 a2i 

B\ 2 a 2 B 12 = e^-^/^-'Wtas, (45) 

where the integer k = or 1 indicates a remaining ambi- 
guity in the phase e tv that we have not yet fixed. This 
implies that the physical loop operators transform as: 

B{ 2 Wx{ai)B 12 = Wi{a x ) 

Bl 2 Wi(h)B 12 = ^■ Kk+iv ^ m - l ''Wl{a 1 )W 1 {b 1 ) (46) 

Using the fact that W^\ai) = TU 2 m ^(6i) = 1, we can 
use the constraint 

[B 12 ,WF-\b 1 )]=Q (47) 

to fix k = 1 in the case when m — I is odd. 

From the above transformation, we can find the braid 
matrix, B\2- First, we pick a basis of the \m — l\ states: 

Wi(6i)|n) = I n + 1 mod \m - l\), 
IUi(ai)|n) =e 2mn ^ m -V\n). (48) 

Then, consider 

B\ 2 \n) = e l6 e l7Tkn+m ' T ' {m - l \wl{a 1 )W l {b l )) n \Q) 

= el e el ^n- l xn 2 /(m-;)| n ) ] ( 49 ) 

where the overall phase e 10 depends on details and is not 
topological. Recall that k = 1 when m — I is odd. When 
771 — Z is even, we cannot fix k using this approach; the two 
different choices of k = or 1 are related by a basis trans- 
formation 1 77) — > 1 77+ (m — I)/ 2). Alternatively, note that 

B 12 (k = 0) = W^ m - l)/2 {b 1 )B 12 {k = l)wt m - l)/2 (h) 
when m — I is even. 

When m — I is odd, we see that the braid matrix B12 
computed in this way agrees precisely with that obtained 
through the bulk geometric approach of Section llll Al (see 
eqn. (fT3j) ). When 777 — I is even, we see that the case 
k = agrees with (fT3[) . 



D. An alternative approach: braiding without 
moving the defects 

Now, instead of carrying out the braiding by a continu- 
ous counterclockwise motion of the twist defects, we will 
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FIG. 9: Illustration of ID protocol for carrying out an ef- 
fective braiding process. This protocol does not require a 
continuous motion in both directions and may be easier to 
implement in many physical realizations, for instance through 
gating. 



consider a somewhat different protocol that can be im- 
plemented by considering a path purely in the ID Hamil- 
tonian of the CFT4£ This protocol may be more readily 
realizable in a physical system by, for example, gating. 
It also helps us compare more directly with the braiding 
statistics of related defects obtained in Ref. [22M25I . 
First, let us define the following Hamiltonian: 



H ab = -\t\(e-^ d ^ m -^ata b 



H.c), (50) 



which couples the zero modes located at the a and b 
defects. This has the effect of generating an energy gap 
in the degenerate subspace formed by the a and b defects. 
The phases 8 a b determine the eigenvalue of ot\oib on the 
ground state. 

Now, we can consider the following two-step process: 



H(r) 



H 2 ^3 = (1 - T)#34 + rff 24 , r € [0, 1] 
ffi_> 4 = (2 - t)H 24 + (t - 1)H 12) t e [1, 2] 

(51) 



As is shown in Fig. [9j in the first half of the process 
t € [0, 1], the Hamiltonian i?2->3 moves the zero mode at 
defect 2 to defect 3. In the second half, the Hamiltonian 
Hi^4 moves the zero mode at defect 1 to 4. We set 
#12 = #34 so that up to a translation, this is a closed 
path in the Hamiltonian projected onto the low energy 
subspace, which exchanges the zero modes ct\ and a 2 . It 
is therefore a one-dimensional protocol for braiding the 
defects that does not require motion in both directions 
and, depending on the realization of the twist defects, 
can be realized physically by gating^ 

In order to understand the effect of these processes 
on the ground state subspace, we first observe that the 



operators 



0\ = a2«3«4, 



2 = a\a 2 a\, (52) 
commute with the two processes ; 22 ' 23 respectively: 



[H2-*, Ox] 
[Hi^a,0 2 ] 



(53) 



Using this, we can now obtain the effect of this process 
on the zero modes ot\ and a 2 . We let V a ^b{T) be the 
projector onto the ground state sector of _ff a ^b(r). First, 
we define integers ki, k 2 , and such that: 

9 34 e {h - 1/2, h + 1/2), 

024 e -{k 2 - 1/2, k 2 + 1/2) - (to - I + l)/2. (54) 
Then, we find 

7W0)Oi7W0) = e 2 ^/^a 2 , 
^3(l)Oi^ 3 (l) = e-" + " (2fe2 - 1)/(m -Vc4a 4 , 



(55) 



while 



7W1)0 2 7W1) 
7W2)0 2 ?W2) 



e -m+in(2k 2 -l)/(m-l) a \ 



1- 



_ i27rfci/(?n-Z) n ,t 



(56) 



Note that a 4 a\a 4 commutes with the second process, 
#2->-3 (t), and so is unchanged. After the process is over, 
we can relabel a 4 as a 2 and 03 as a±. Therefore, we 
conclude 



B\ 2 a 2 B 12 = -e l 
B\ 2ai B 12 -- 



r (2fc 2 -2fci-l) t 

• y >a 2 a[a 2 , 



a i-n_ T (2/s 2 -2fei-l) 



a 2 . 



(57) 



From this, and following the same steps as in the previous 
section, we find that the ID braiding protocol gives 



B 12 \n) 



e i6 e i-K{n-(k 2 ~k 1 )) /(m-Z)| n \ 



(58) 



Comparing with (14U)) , we see that this ID braiding pro- 
tocol yields an extra integer degree of freedom, [k 2 — fci), 
which depends on the phases 9 24 and 6*34. We can in- 
terpret these braidings as the pure braiding of (|4T)1) . 
combined with the action of the Wilson loop operators 



IV. Zi CHARGE CONJUGATION TWIST 
DEFECTS IN SINGLE-COMPONENT ABELIAN 
STATES 

A. Charge conjugation twist defects 



Ref. I22M25I considered another type of extrinsic de- 
fect created by the coupling between a superconductor 
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and the edge states of a fractional quantum Hall state 
or a fractional quantum spin Hall state. In a fractional 
quantum Hall state, the defect occurs between the nor- 
mal inter-edge tunneling region and the superconducting 
pairing region. In a fractional quantum spin Hall state, 
the defect occurs between superconducting and ferromag- 
netic regions along the edge. In both systems, the role of 
the superconductor is to break charge conservation. In 
the following we will not review their proposals in more 
detail, but we will discuss the understanding of this de- 
fect as a twist defect, and its relation to the Z 2 twist 
defect discussed in the previous section. 

To begin with, we consider Abelian topological states 
that are described by U(1)n CS theory, namely the 1/N- 
Laughlin FQH states: 



N x 
C = — e ^%cW 
Air 



(59) 



These states have N quasiparticles, {7.;} for i = 1, • • • , N, 
with topological spin 9 a and fractional charge q a 



ira 



q a = a/N. 



(60) 



9 a is defined modulo ir when TV is odd (fermions) and 
modulo 2tt when N is even (bosons), while q a is defined 
modulo 1. The mutual statistics between quasiparticles 
7 a and j b is given by 



lb = 2irab/N, 



(61) 



which is defined modulo 2ir. 

We see that when charge conservation is broken, this 
theory has a Z2 symmetry associated with 

la ^ IN -a- (62) 

In the CS field theory, this symmetry is implemented as 

a — > —a. (63) 

This implies that we can consider twist defects associ- 
ated with this Zi symmetry. As in the previous section, 
we align all the twist defects, and we cut the system 
along a line to obtain counterpropagating edge states. 
The electron operator on each edge is: 



cL 



J N 4>L 



eR 



-iN<j> R 



(64) 



Now, in the untwisted A regions, we consider the usual 
hopping, ^l L ^ e R + H.c. oc cos(A0). In the twisted 
regions, we apply the Z2 action to one of the chiral 
edge states: 4>r — > — </>_r, to get the tunneling term 
^Il^Ir + H c - x cos(jV0), where 



5R, 



(65) 



Thus we have: 



6H t = Z 



(^i + H.c.) if xeA t 
(4^ + tf.c.) if xeB l 



cos(N(j>) 
cos{N6) 



if x e Ai 
if x e Bi 



(66) 



This is precisely what has been considered in Ref. I22M24 
[27l . where it was observed that the defects have quantum 
dimension ^/2N. 



B. Relation to genons in two-component states 

Now let us consider the Z2 genons in the two- 

m I 



component theories with .ftT-matrix K 



I 



In 



this theory, the two kinds of tunneling terms that we 
add are (see eqn. (J2SJ)) 



SXJ I 171 + 1 A " 

otit — g cos — - — cp + 



cos(^0_) if x&Ai 
cos(2^i^_) if xeBi 



Now observe that if we define 



*>R± 



b L l±< 

i>Ri ± 1 



•L2, 



and further define 



then 



6W.t = 9 cos 



(67) 



(68) 



(69) 




if x € A t 
if x e Bi 



(70) 

Since both tunneling terms contain cos( ii y^0+), <fi+ is 
pinned everywhere and in the low energy sector can be 
replaced by a constant. Now we see that in the Luttinger 
liquid theory, the two different tunneling terms are iden- 
tical to the superconductor - ferromagnetic terms in the 
FQSH edge (see (|6"6"|) ). This explains the agreement be- 
tween the braiding matrices computed from the Dchn 
twists of the high genus surface ( Section IIII[) . and that 
computed in Ref. [22H241 

In (|66p , the defects have quantum dimension d = \/2N, 
while in the two-component case, they have quantum di- 
mension d — \J\m — l\. Therefore, d 2 is not restricted to 
be even in the latter case. The reason is that in the two- 
component case, the coefficient of the boson fields in the 
cos terms can have half-integer values (see [70]) . while it is 
restricted to be integer in the FQSH setup (1551) . This can 
be traced to the fact that the boson fields (9_ and 0_ are 
compactified on a circle: (6L, </>_) ~ (#_ + 2tt, 0_ + 2ir). 
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Therefore, cos( ii y^6'_) and cos( I!2 j^0_) are not individ- 
ually invariant under such a gauge transformation if m — l 
is odd. However, in the two-component case, this is al- 
lowed, because whenever one of 4>li or (f>m are advanced 
by 27r, both and 0+ or 0_ and 0+ will change by 
2ir, so that the products cos( I2 j^</> + ) cos(^^O-) and 
cos( !i Y^(/) + ) cos(2^0_) are invariant. 

In light of this, we note that even in the n-layer FQH 
states, one can introduce superconductivity in addition 
to twisted tunnelings in order to create zero modes. We 
expect then that one can always map such situations onto 
a 2n-layer system without superconductivity but with 
twisted tunnelings. 

V. Z 2 GENONS IN TWO-COMPONENT 
NON-ABELIAN STATES AND UNIVERSAL 
QUANTUM COMPUTING 



7r/8 phase gate: 



(2) , a m , a 




A large part of the analysis in the last section can 
be generalized to non-Abelian states. In this section we 
briefly comment on the somewhat more exotic possibil- 
ity of twist defects in non-abelian states. Let G denote 
any non-abelian topological state, and let us consider two 
independent copies of such a state, which we will label 
G x G. For example, we may take G to be an Ising 
topologically ordered state. The G x G state has a Z 2 
symmetry in its topological quantum numbers associated 
with exchanging the two copies, n pairs of genons on a 
sphere will therefore lead to a single copy of G, on a 
genus g = n — 1 surface. The ground state degeneracy 
S g is then given by the general formula^ 9 -^ 



S a 



D 



N„ 
i=l 



2(9-1) 



(71) 



where di is the quantum dimension of the ith quasi- 
particle, N qp is the number of quasiparticles in G, and 



D = \JJ2i^i d% i s tne total quantum dimension of G. In 
the limit of large g, we see that S g ~ D 2n , which shows 
that the genons have a quantum dimension 



dz, = D. 



(72) 



In the case of genons in the Ising x Ising theory, we 
see that dz 2 = 2, and can effectively yield a single Ising 
theory on a high genus surface. The braiding of the twist 
defects corresponds to Dehn twists of the Ising theory on 
the high genus surface. 

While the Ising theory by itself is known to be non- 
universal for TQC, it is known that it can be made uni- 
versal if it were possible to add handles in the space and 
carry out Dehn twist a 29 i 30 . The Z 2 genons provide a 
physical realization of such a possibility and therefore 
can be used to render universal for TQC a state that 
without the twist defects is non-universal for TQC. 

Here we will briefly describe how to use the genons 
to implement the necessary computational gates in the 



FIG. 10: Illustration of operations (1), (2), (3), and (3'), de- 
scribed in the main text for performing the CNOT and 7r/8 
phase gates. The red circles indicate that the topological 
charge of the two twist defects enclosed in each circle is trivial. 
The black loops with arrows stand for quasi-particle paths. 



Ising x Ising state, adapting the basic ideas of Ref. 
I29II30I . For more comprehensive reviews of the Ising topo- 
logical state, we refer the readers to books and review 
article o 2 ' 50 . The Ising x Ising state contains 9 distinct 
quasiparticles, of the form a x b, where a, b = l,ip,<j are 
the three quasiparticles of the Ising state, which have 
spins hi = 0, hi/, = 1/2 and h a = 1/16. In the pres- 
ence of genons, the non-contractible loops surrounding 
the genons can be mapped to non-contractible loops of 
a single Ising theory on a high genus surface; therefore 
each non-contractible loop surrounding the genons can 
have one of three topological charges: 1, ip, or a. 

A universal set of gates for quantum computing is ob- 
tained in terms of the single-qubit '71-/8' phase gate; 29 i 30 



Gi = 

and two two-qubit gates: 

/l \ 
10 
10 
\o -1/ 



1 







e 2 ™/ 8 



(73) 



G 2 



G, = 
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\r* 








1/ 



(74) 



G2 is, after a basis transformation, the controlled-not 
(CNOT) gatei 3 ^ Each pair of a quasiparticles in the Ising 
state can fuse to either 1 or ip, and thus these two possible 
states form the two states of a qubit. Therefore n qubits 
requires In a particles. 

While G3 can be implemented by braiding the quasi- 
particles in the Ising theory, Gi and G2 cannot be implc- 
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mented by braiding quasiparticles alone. However they 
can be implemented by creating genons and braiding 
them around the quasiparticles. In order to implement 
G\, we consider two quasiparticles of type a x 1, which 
we label 1 and 2. The pair (1,2) can be either in the 
lxl channel, or the -0x1 channel. We implement the 
following processes (see Fig. [TU|) : 

(1) Create two pairs of genons out of the vacuum, 
which we label gi,...,^. Let C gu9i+1 label the loop 
surrounding gi and <?;+i. After creating the genons, en- 
sure that the topological charges around the loops C gi „ 2 
and C ff3jff4 are trivial, which can be done by performing 
a measurement. 

(2) Braid quasiparticle 2 around C g2t93 . 

(3) Perform a double-exchange of g 2 and (73. 

(4) Re-annihilate the genons. 

This procedure implements G\. After step (2), the 
topological charge of C gil g 2 is a. It can be shown that 
if (1,2) are in the lxl channel, then the topological 
charge of C g2t93 is (1 — ^)/v / 2 - 29 ' 30 Therefore the double 
exchange, g 2 and 53, which implements the double Dehn 
twist in the geometrical picture, has an eigenvalue of 1, 
because e 2 ^ 7 "' 11 ' = e 2{2lTlh ^ = 1 . On the other hand, if 
(1,2) are in the -0x1 sector, then the topological charge 
of C g2i g 3 is a £2- In this case, the double exchange of the 
genons gives a phase corresponding to twice the spin of 

a . (g27ri/16)2 _ e 27rj/8^ 

In order to implement G 2 , consider two pairs of quasi- 
particles, all of type ax 1, which we label 1, 2, 3, 4. Then, 
perform the same steps (1), (2), and (4) as above replac- 
ing (3) with (3'): 

(3') Braid the pair (3,4) around the loop C g2 , g3 . 

Based on the previous discussion, the step (3') gives a 
+1, unless (3, 4) is in the state ip X 1 and the topological 
charge of C g2 _ g3 is a (ie if (1, 2) is also in the state ip x 1), 
in which case step (3') gives —1. 



VI. Z z TWIST DEFECTS AND BRAIDING IN 
THREE-COMPONENT ABELIAN FQH STATES 

To understand more generic behavior of twist defects, 
here we extend the analysis of quantum dimensions and 
braiding of twist defects to Z% twist defects in a three- 
component Abelian state, which are described by Z7(l) x 
[7(1) x U(l) CS theory: 



C 



4tt 



(75) 



for a rank-3 ii"-matrix. Such Z3 defects can be realized by 
the lattice dislocations in the Chern number 3 topolog- 
ical nematic states*. In principle it can also be realized 
in triple layer quantum Hall states, which have also been 
realized experimentall y 51 ' 52 . Imposing a Z3 layer symme- 
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tr 



T 



T4- 



— !— 



FIG. 11: Oriented all of the defects along a line and cutting 
yields three counterpropagating edge modes. Since the defects 
are Z3, there are now three different regions, A iy Bi, and d, 
where we introduce different twisted tunneling terms. 



try implies that the A'-matrix depends on two integers: 



K = 



I 



* m I 

I m I 
\ I I m 1 



(76) 



As we will see, this case requires a more general anal- 
ysis with features that do not show up in the case of 
the simpler Z 2 twist defects. Since the bulk geometric 
picture and mapping to higher-genus surface is harder 
to visualize than the Z 2 case, we will first carry out our 
analysis based on the 1+1D edge CFT picture, and then 
present the bulk geometric picture. 



A. 1 + 1D edge CFT picture of the Z 3 defects 

In this case, the Z$ twist defects come in groups of 
three. As before, we align them on a torus and cut the 
system along the line to obtain counterpropagating edge 
states. Similar to the Z2 case, the actions of the edge the- 
ory on the two sides of the cut are given by ([2"Tj) and l[2"7j) 
respectively. As shown before, we can obtain the commu- 
tator ((24)) and (|28|) by quantizing the theory. However, 
instead of using (|3"Tj) , here we choose a different scheme 
to take care of the correct statistics between electrons on 
different layers. We write the electron operators as: 



^eLI = VLie 



eRI 



iKu4>LZ 
-iKu<pRi 



(77) 



where t]li and rjm are Klein factors that satisfy 

Vli = Vri = 1, Vli = Vli, Vri = Vm, 

I [vli, T)lj} = [vri, VR.A = [vli,Vrj] = 0, for m-le 
\ {vli,Vlj} = {vri, VrA = {Vli, Vrj} = 0, for m - I o 

(78) 

When we glue the edges back together, we now have three 
different regions where we introduce different tunneling 
operators (see Fig. ITTj) : 



8H t 



tun 




K-i cos(Kij(j)j) if x € Ai 
K.' I cos(Kij(f)'j) if x e Bi 
K'lcasiKijtf'j) if xeG 



(79) 
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where now 



<t>'l = 4>LI + 0R(J+1)%3) <Pl = 4>LI + <t>R(I+2)%3, 



■m — l 



■ra — l 



VLIVR(I+1)%3i 



1C" = i m l VLlVR.{I+2)%3- 



(80) 



The factor i m ~ l helps keep the hermiticity of the opera- 
tors K,i, IC'j and K![. In the absence of the twist defects, 
we only consider a single A region, and the minima of 
the cosine potential give |Det K\ = |(m + 2l)(m — ^) 2 | 
different states. These states can be labelled by integer 
vectors p~ (p+,Pi,—2,i,Pi,o,—i), such that 



0*^(1,-2,1; 



\P) =e^w|p), 

' 27r P(l,-2,l) 

\p) = e ™-> \p), 



e^d,o,-i) |p) = e s=i 



(81) 



where we have defined 



a(f>i + bfo + efts, and 



(a,b,c) 

0+ = 0(i,i,i)- 

Now consider the system with n > triplets of twist 
defects. While formally the tunneling Hamiltonian in 
the different regions differ only by a cyclic permutation 
of the layers, there is an important difference between the 
regions. Since all physical operators can only be built out 
of electron operators, there is a gauge symmetry in the 
values of 4>li and 4>ri- The following transformation: 



Hi -> 

f>RI -> 



•>L1 



- 2-KKjjrij, 



(82) 



with m, n2, «3 £ Z, preserves all physical opera- 
tors and are thus considered gauge symmetries of the 
theory. All physical operators must be gauge invari- 
ant. In the A regions, the quasiparticle tunneling op- 
erators e 1 ^ 1 ^^' are all gauge- invariant. In contrast, 
in the B and C regions, e 1 ^ 1 and e 1 ^ 1 are unphys- 

ical. However, the operators e ^+^ B >' = e l ^+^ XB i\ 
e #"(x 0j ) _ e i<t>+{x Ci ) ^ e #'(i,o,-i)( a: Bi)+^'(-i,i,o)( ;,: c i ) anc j 

e ^(-i,i,o)(^B i )+0'('o,-i,i)(*c i ) are alg0 physical, gauge- 
invariant operators. The cosine potentials pin the eigen- 
values of these operators. Due to the commutation rela- 
tions: 

[(/>i(x), 4>'j{y)\ = i7T ( K 7.lj+i)%3 ~ K 7,j) a s^( x - y)> 
[<f>i(x), 4>j{y)] = i7r ( K 7}j +2)%3 - K i}) s s n ( x - v): 

Wj(x),^(y)] = in(K-l +1)%3 - Ky}) S ga{x - y), (83) 

we see that the largest set of independent commuting 
physical operators is e i<t,l{x ^\ e ^+<> B i), e #+(*c i ) j anc i 
e i <f>{i,o,-i)i x B i )+i<P['-i,i,o)(xa i ) 

Minimizing the cosine potential in the A regions, we 
get |(m + 2l)(m — l) 2 \ states for each A region, as in the 
case with no twist defects. In the B and C regions, the 
operator e l< ^+ is pinned to one of (m + 21) values in each 
region. Finally, for each triplet of defects, the operator 



e^'d.o.-D^B^+^'-i.i.o)^,) is pinne d to onc of | m _j| va i_ 

ues. In Appendix [AJ we will provide a more detailed 
treatment which explicitly includes the Klein factors. 

The above counting gives a total of \(m+2l) 3n (m— l) 3n \ 
states. However, there are additional constraints: 9 The 
electric charge is a local observable, and therefore states 
with different charges at the defects are not topologically 
degenerate. This gives 3n — 1 constraints on each of the 
eigenvalues of e"^ + in the different regions. The total 
number of states is therefore \(m + 2l)(m — l) 3n \. Each 
Z3 twist defect therefore has quantum dimension 



dz, = \m-l\ 



(84) 



Writing the ground state degeneracy as |Det K\\m — 
l\ 3n ~ 2 , the factor |Det K\ can be understood from the 
fact that we started on a torus, while the factor \m — 
l\ 3n ~ 2 can be understood as due to the appearance of a 
non-trivial quasiparticle loop algebra, which forms 3n — 2 
copies of a magnetic algebra. The degeneracy \m — l\ 3n ~ 2 
is equivalent to what would be obtained in a U(l) m -i 
CS theory on a genus g = 'in — 2 surface (the explicit 
mapping to the high genus surface will be explained in 
Sec. IVI Bp . Below, we calculate the braiding of these 
Z% genons, which can also be related to Dehn twists of 
[/(l) m _z CS theory. 

To compute the braiding, as before, we can construct 
zero modes localized to each domain wall, which corre- 
spond to quasiparticle tunneling around the defects: 



mRV2Re 

r\lRT\2RZ 



Vir^rs 



i0i(xA i ) e -#' 1 (x Bj ) ) i m od 3 = 1 
i*3(^ 4 ) e -»y 8 '(xo 4 ) i i mod 3 = 2 

= 



i^(x Ci ) e -iMxA t+1 ) ■ mod 3 



(85) 



Similarly, we define : 

ri2RV3Re i4 ' 2ixA ^e" l ' t ''^ XB ^ ) , i mod 3 

Pi={ mRmR^'^^e-^"^, i mod 3 

V2RV3Re i ^ (xc ' ) e~ l ' t ' :>{xAt+1 \ i mod 3 



(86) 



and similarly for 7,;. Fig. [T2l displays the quasiparticle 
tunneling process described by these operators. For ex- 
ample, cti describes the process in which a layer 2 quasi- 
particle tunnels around a twist defect and goes to layer 
1. The Klein factors are added so that these operators 
are zero modes: 

[8H tm ,oLi] = [SH tun ,/3i] = [SH tun ,li\ = 0. (87) 

Since = 4>\ + cf>2 + 4>3 is fixed everywhere, there is a 
local constraint between the three zero modes similar to 
that in the Z2 case (c/. eq. (l38|) ): 



a k (3 k jk = e—'e lQk , 



where Qk = ^((j)+(xk) — (f>+(xk-i)) is the charge on the 
kth defect, where here Xk refers to the region between the 
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k — 1th and kth. defect. In what follows, for simplicity 
we set Qk = 0, as this does not affect the topological 
properties of the defects. 



a n a n+ k = a n +kOin&"^ 1 ', Pnfin+k = 

O-nPn+k = (-l) m_i /3n+fc a ne^ rrT , finln+k = 

a n 7n+k = (-l)" l ~ i 7n+fc Q! ne" I ^ T , Pn a n+k = 

a n (3 n = (-lr-'^^e^, /3„ 7 n = (- 
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3i+1 3i-t-2 3i+3 

II 

a 3i,l f X-| - - - - \-X I f x I *3i + 3 

FIG. 12: Physical meaning of the zero mode operators at, pi 
and 7i for the Z3 twist defect case. In subsequent figures, we 
will use the bottom figure, which is more precisely described 
by the top figure. 



where k > 0. The first line shows the algebra between 
the same zero mode operators located on different de- 
fects. The algebra of different operators on different sites 
is shown on the second and the third line, while in the 
last line, we write down the algebra between different op- 
erators on the same defect. The zero modes can be used 
to construct the quasiparticle Wilson loop operators: 

W(a 3 k+l) = a 3 k+lf3 3 k+ 2 J 3 k+3, 
W(b 3k +l) = 73/c+3a3fe+4/33fe+5, 
W(a 3 k+ 2 ) = /33k+lj3k+2&3k+3, 

W(b 3k+2 ) = a 3+ 3/3 3fe+ 473/c+5, 

W{a 3k+3 ) = ~e^a\ k+2 a 3k+3 , 

W(b 3k+3 ) = -e^ pl k+2 f3 3k+3 , (90) 

where we have chosen the overall phases so that (see Ap- 
pendix [B] ) 

W{a t ) rn - 1 = W{bi) m ~ l = 1. (91) 

Fig. [T3] displays the loops a* and bi along which the 
quasiparticle tunnels by the action of the Wilson loop 



Using the definition of the zero modes, we obtain the 
following algebra: 



r 




FIG. 13: Physical meaning of the quasiparticle Wilson loop 
operators. 

operators W(a,i) or W(bi). These operators satisfy: 

W{a 3k+l )W(b 3k+j ) = W(b 3k+J )W(a 3k+l )e s » 2 ™/( m ~ l \ 

(92) 

which leads to \m — l\ 3n states forming the irreducible 
representation of this algebra. Adding an extra factor of 
\m + 2l\ for the possible values of e i(<#,1+<fc+<fe) , we get 
a total of |(m + 21) (m — l) 3n \ topologically degeneracy 
states. This confirms again that the quantum dimension 
of each Z 3 defect is dz 3 = \m — l\. 

In order to calculate the braiding matrices, we need to 
understand how the zero mode operators transform under 
a braiding process. From Fig. [TJJ we see that under a 
clockwise braiding of 2 and 3, B2 3 a2B 23 = e l ^a 3 and 
B 23 a 3 Bl 3 = e l6 jl 720:3, where B 23 is the braid matrix. 
The Z 3 layer cyclic symmetry implies that B 23 j3 2 B\ 3 = 
e^/? 3 , B 23 (3 3 Bl 3 = e l6 a\a 2 p 3 , B 23l2 B\ z = e^ 73 and 
lz B\ 3 = e l9 pl(3 2 ~/ 3 . We can fix = e~ ie by using the 
fact that the braiding operation should commute with 
the loop operator W{a\) and W{a 2 ) that encircle these 
two defects. Using (1551) . we then have = 1. Thus, 



I 

Pn+fePne™-' , 7„7n+fc = "fn+kj n e m -' , 

(-l) m ~ l "f n +kl3 n e~ , j n a n+ k = (-l) m ~ l ~fn+ka n e^ , 
(-l) m - l p n+k a n e^, ln p n+k = (-l) m - l (3 n+ kl n e^, 
l) m - l ln(3 n e^ , j n a n = (-l) m - l a n7n e i ^. (89) 
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FIG. 14: Effect of a clockwise braid on the zero modes on the second and third dislocation, (a) Q2 — ► «3 (b) «3 gets transformed 
into the combination of 7372 and 03. Therefore, as —¥ 737203. Similar to the Z2 case, the ordering ambiguity or additional 
phase factor can be fixed with additional constraints. 





FIG. 15: Illustration of the ID protocol for carrying out an effective braiding process of the Z3 twist defects. The purple, 
orange and gray circles mark up the zero modes located on the corresponding defects. 



we find: 

B 23 a 2 Bl 3 = e^O!3, B 23 a 3 Bl 3 = 7372^3, 
B 23 /3 2 b1 3 = e^/3 3 , B 23 (3 3 Bl 3 = e^^ala 2 /3 3) 

^2372^3 = e^ 73 , ,82373^3 = Plfalz, 

(93) 

where the integer k = 0, 1, 2 is a phase ambiguity that is 
not yet fixed. Further, considering (l91"j) , we obtain k = 
when m — I is not a multiple of 3. For m — Z being a 
multiple of 3, we still have 3 different choices of k which 
is similar to the Z 2 case. 

Now we choose a basis for the \m — l\ 3n different states, 
|{rii}), for rii = 1, • • • , |m — l\ so that 

W(a i )|{n J -}) = e 2 ' ri "*/( ro -')|{n J -}) ) 
W r (6 i )|{» i }) = |{(n i +i ii )%|m-l|}). (94) 

Note that we can ignore the extra degeneracy of \m + 2l\ 
associated with the possible eigenvalues of e l ^+ , as they 



are independent of the braiding. In this basis, we find 
that the braid matrix is: 

({n^\B 23 \{ ni }) 

The braid matrix B 23 can be viewed as an element of 
the mapping class group in the Z7(l) m _; CS theory on 
a high genus surface. Observe that -B23 only has a non- 
trivial action on an \m — ^-dimensional subspace of the 
ground states, associated with the states that form the 
\m — ^|-dimesional irreducible representation of W(a 3 ) 
and W(b 3 ). Now consider the Dehn twists U a and Ub, 
around the a and b cycles of the torus for J7(l) m _/ CS 
theory (see JTSJ), dHJ), (fT7])V When m - I is odd, the 
action of -B23 in this relevant subspace is 

B 23 = e w {U a U b )\ (96) 
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where the phase e depends on details of the path. When 
m — Z is even, 

B 23 = e w W(b^ m - l)/ \U a U b )\w\bz)) (m - l)/2 (97) 

Therefore, within this subspace, B 23 coincides with a se- 
quence of Dehn twists in J7(l) m _j CS. When f? 2 3 is not 
projected onto this subspace, it can be viewed as a se- 
quence of Dehn twists in U{l) m -i OS theory on a high 
genus surface. We will also see this result via the geo- 
metric construction of the subsequent section. 

Similar to the Z 2 case, we can also realize the braid- 
ing through a purely ID protocol, without continuously 
moving the defects in both dimensions. First, we define 
the following Hamiltonian: 



Habc = -\t\(e- 



m - 1 a a p b "/ c + e 



"aba it + 

=i e.m-i a a Oib + H.c.j, 



which couples the zero modes located at the a, b and c 
defects. This Hamiltonian will gap out the degenerate 
subspace formed by the a, b and c defects. 
Now, we can consider the following processes: 



H(t) 



H 3 ^ 5 = (1 - t)H 456 + tH 346 , t e [0, 1] 
H 2 ^6 = (2 - r)H 3i6 + (r - 1)H 234 , r G [1, 2] 

(99) 

As is shown in Fig. [15l in the first half of the process 
t e [0,1], the Hamiltonian H 3 ^ takes the zero modes 
on defect 3 to defect 5. In the second half re [1,2], the 
Hamiltonian H 2 ^q takes the zero modes on defect 2 to 
6. We set #456 = #234 so that up to a translation, this 
is a closed path in the Hamiltonian projected to the low 
energy subspace. This process effectively exchanges the 
zero modes at defects 2 and 3. 

In order to understand the effect of these processes on 
the ground state subspace, we first observe that operators 

Ox = a 3 /3 4 ~f 5 a 6 , 0[ = J 3 a 4 f35l6, 

02 = a 2 /3 3 j 4 a 6 , 0' 2 = 720(3,0476 (100) 

commute with the two processes, respectively: 

[fls-ns.Oi] = [H^,0' x ] =0 

[H 2 ^ 6 ,0 2 ] = [H 2 ^ 6 ,0' 2 ] =0 (101) 

Using this, we can obtain the effect of these processes on 
the zero modes. We let V a ^b(T) be the projector onto 
the ground state sector of H a ^b. First, we define integers 
k\ and k 2 such that 

0456 € (fci-l/2,fci + l/2), 
#346 e (fc 2 - 1/2, fc 2 + 1/2). (102) 
Then, we find that 



TVus (OjdTV* (0) = e 
^5(0)01^5(0) 
7W1W 3 - 5 (1) =e 
7^5(1)0^5(1) =e 



-' a 3 , 



e ™-' 73, 

i27Tfc, 4. 



m_! 7e75Q:6, 
~^4&76, 



where we have used the relations (a 3 (34je) m 1 = 

(/3 3 7406) m ~' = 1 (see Appendix |B|) and OiA'7i = e m ~ l 
(assuming the charge on each dislocation is 0). Also, we 
have 

V 2 ^e(l)0 2 V 2 ^ 6 (l) = 
^6(1)02^6(1) = e^ 72 , 
^6(2)02^6(2) = e^a 6 , 
^6(2)02^6(2) = e^ i 76 . (104) 

2iii 

Therefore, with the relations oiiPiji — e m -' , we conclude 

. i27r(fc 1 -fc2) 

B 23 a 2 B 23 = e a 3, 

, i2ir(fc 2 -fcl) 4. 

B 23 a 3 B 23 = e ™-' 73720:3, 

+ i27r(fci + 2fc 2 ) 

B 23 f3 2 Bl 3 = e —< /3 3 , 

, i2-7,(2k 1 +k2) . 

B 23 /3 3 B 23 = e ™-i a 3 a 2 /3 3 , 
B 23 l 2 B\ 3 = e m-i 



73, 



£2373^3 = e" 2 /3|/3 273 . (105) 

This result is slightly different from (|93l) . The reason is 
that when we restrict the discussion to the same charge 
sector, in which the relations aiftiji — e m -' hold for any 
i (see (|88p) , the Hamiltonian (|99p of this ID protocol 
breaks the Z 3 cyclic layer symmetry for generic k± and 
k 2 . If m — I is not a multiple of 3, only when k-y = k 2 = 
is the Z 3 layer symmetry restored. If m — I is a mul- 
tiple of 3, the condition of the symmetry is k\ = and 
k 2 = 0, ±(m — 0/3- Once the symmetry is restored. This 
result agrees with (l93|) . On the other hand, if we choose 
to loosen the charge constraints but insist on the Z 3 layer 
symmetry of (|99l) (by adding other terms related to the 
original terms by cyclic symmetry into the definition of 
(|98p). then the charge on each dislocation changes gener- 
ically during the "braiding" process dictated by the ID 
protocol Hamiltonian (|99|) . 



(103) 



B. Bulk geometric picture 

When I = in the K matrix in Eq. ([75| . the three 
layers are decoupled and it is possible to understand the 
Z 3 twist defects by directly mapping the tri-layer system 
with twist defects to a high genus surface, similar to the 
bilayer Z 2 case. The mapping is probably too compli- 
cated for practical purposes, and we will only discuss it 
in the simplest case of a sphere with three Z 3 defects, 
as an illustration of the general situation (as in the Z 2 
case, we expect the mapping to also be possible when 
I 7^ 0, though this is even more complicated). Illustra- 
tion of the mapping is given in Fig. [T6j We consider a 
cut-and-glue procedure to describe the defects. Cutting 
the sphere along a line connecting the defects 1,2,3 leads 
to chiral edge states of the three layers around the cut, as 
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FIG. 16: Illustration of the equivalence between a sphere with three Z3 twist defects and a torus, (a) and (b) shows the sphere 
with twist defects cut along the branch-cut line connecting defects 1, 2, 3. The red, blue and green lines denote the edge states 
of the three layers along the cut. After the cutting, the system in (b) is equivalent to three decoupled systems of the three 
layers shown in (c). Then the three layers are glued along the branch-cut line after a twist. Along line 1 — 2 the three layers 
are glued in the pattern of red-to-blue (rb), blue-to-green (bg) and green-to-red (gr), in which the first and second colors are 
the layer index of the lower and upper edges respectively. The result of this gluing is shown in (d) which is equivalent to (e). 
Then the green and blue layers are glued along 2 — 3 line, shown in panels (f) and (g). Finally, the rest of the edge states in 
(g) are glued following the rule of red-to-green, green-to-blue, which leads to the torus shown in (h) and (i). 



is shown in Fig. [16] (a) where the three layers are labeled 
by red, blue and green lines. As is illustrated in Fig. Ql)] 
(b) , the system with the cut is equivalent to three decou- 
pled layers, each of which is topologically a disk. Then 
the twist defects are created by gluing different layers 
shifted by a Z3 operation. Along the branch-cut between 
defects 1 and 2 the left-movers on the lower edge in the 




FIG. 17: A loop around two Z3 twist defects 2 and 3 (the 
dotted lines) on the sphere (a) is topologically equivalent to 
a nontrivial loop on the torus shown in (c). This is the same 
loop that appears in Fig. [14] (b). The color of the dotted 
line indicates the layer indices. This equivalence is obtained 
by following the steps shown in Fig. 1161 Panel (b) shows 
an intermediate steps of the topological deformation which 
corresponds to Fig. [16] (d). Other intermediate steps are 
omitted. 



red (r), blue (b) and green (g) layers are glued with the 
right movers on the upper edge in the b,g,r layers, re- 
spectively. Along the branch-cut between 2 and 3, the 
layers are glued in the opposite fashion so that r,b,g left- 
movers are glued to g,r,b right-movers, respectively. This 
gluing can be done in three steps. First, the layers are 
glued along the branch-cut between 1 and 2, which glues 
the three disks into one disk, as is shown in Fig. [16] (d) 
and (e). Secondly, along the branch-cut between 2 and 3 
the green left-mover was glued to the blue right-mover, 
which leads to a sphere with two punctures shown in Fig. 
1161 (g) . Thirdly, the remaining red and blue left-moving 
edge states are glued to green and red right-movers, re- 
spectively, which glues the punctured sphere to a torus. 



We can also see from the geometric picture that the 
Wilson loop operators indeed correspond to a non-trivial 
loops on the torus. For example, the loop correspond- 
ing to the operator 7^73 is shown in Fig. [T7] (a). In 
the mapping this loop is mapped to a nontrivial loop on 
the torus, as is illustrated in Fig. [T7] (b) and (c). The 
braiding of the defect 2 and 3 will deform this loop to an- 
other non-trivial loop on the torus. Thus, we know that 
the braiding of the defects correspond to a non-trivial 
element of the mapping class group. 
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VII. DECONFINING THE GENONS: 
ORBIFOLD STATES 

So far, we have considered twist defects (genons), la- 
beled by a symmetry group element g £ G, where G 
is a symmetry of the topological quantum numbers of a 
topological phase. These genons are extrinsic defects of 
the system, and separating them costs energy that grows 
with the distance between them, similar to vortices in 
a superfluid. It is possible to deconfine these defects to 
turn them into intrinsic quasiparticles of a neighboring 
topological phase. One way to do this is by finding a way 
to gauge the symmetry G. Such a gauging process corre- 
sponds to orbifolding in the edge CFT^ In this section 
we will discuss the relation between the braid matrices 
of the genons in the confined (ungauged) and deconfined 
(gauged) states. 

In the case of Z 2 twist defects studied in Section Hill 
one can gauge the Z 2 symmetry associated with ex- 
changing the two [7(1 ) gauge fields. Physically, we ex- 
pect that this can be done either by proliferating double 
twist defects^ or by condensing anyons 54 . This leads 
to the U(l) x 17(1) x. Z 2 CS theory studied in Ref. H, 
which describe a set of novel orbifold FQH states^ 2 - and 
"twisted" Z n topological states^ There, the Z 2 twist 
defects simply become finite-energy Z 2 vortices, which 
carry quantum dimension \/\m — l\. These states con- 
tain an Abelian quasiparticle that carries the Z 2 gauge 
charge; when it condenses, the system undergoes a Z 2 
Higgs transition in the 3D Ising universality class, and 
we obtain the Abelian U(l) x U(l) theories. 

There is a close relation between the braiding of the 
Z 2 vortices in the orbifold states and the braiding of the 
extrinsic defects in the Abelian states. To understand 
this, let us briefly review the properties of the orbifold 
states. We recall that in the U(l) x U(l) CS theory, 
for n pairs of Z 2 twist defects on a sphere, the ground 
state degeneracy is |m — ZI™" 1 . When the Z 2 symmetry 
is gauged, a n of these states are Z 2 invariant, while /?„ 
are Z 2 non-invariant, where a„ + f3 n = \m — l\ n ~ , and 
it was found 26 that 



I3n = 



(|m-Z|"- 1 +2"- 1 )/2 
(|m-/|"- 1 + l)/2 



\m- /I"" 1 - 2"- 1 )/2 
dm-/!"" 1 - l)/2 



for 
for 



for 
for 



\m - 
|m 



l\ even 
-l\ odd 



(106) 



\m — l\ even 
|m — l\ odd 



(107) 



Therefore, in the orbifold/twisted states, where we 
gauge the Z 2 , the number of ground states in the pres- 
ence of n pairs of Z 2 vortices is a n . This result allows us 
to deduce the fusion rule o 26 i 32 



( 7 x 7) B = a n l + f3 n (/> + 



(108) 



where <j) is the Abelian quasiparticle that carries the Z 2 
gauge charge, 7 represents the Z 2 vortex, and 7 its con- 
jugate. Therefore, the braid matrix for the Z 2 vortices 
in this case will have dimension a n instead of |m — Z| n . 
Thus for \m — l\ > 2, the braid matrices for genons in 
the confined and deconfined states are not the same, and 
even have a different dimension. 

As we found in Section IIII1 the braiding of the genons 
corresponds to Dehn twists of the U(l) m -i CS theory. 
Under the action of the Dehn twists, Z 2 invariant states 
are mapped to Z 2 invariant states, and similarly for the 
Z 2 non-invariant states. Therefore, the \m — Z|" _1 di- 
mensional braid matrix of the genons, when restricted 
to the Z 2 invariant subspace, will be equivalent to the 
a„-dimensional braid matrix of the Z 2 vortices. 

In the single-component case, described by U(1)n CS 
theory, we considered Z 2 twist defects associated with 
a — > —a, which give defects with quantum dimension 
%/2A. We can consider gauging the Z 2 symmetry in this 
theory, to obtain a CS theory with disconnected gauge 
group U(l) x Z 2 = 0(2)^. The edge CFT of such a 
CS theory is the U{1) N /Z 2 orbifold theory^s When 
N is even, the Z 2 vortices in this theory have quantum 
dimension y/N/2, while when N is odd, they have quan- 
tum dimension yN. This provides an example where 
gauging the symmetry associated with the twist defect 
appears to also change its quantum dimension! 

In the Z 3 case, gauging the Z 3 symmetry leads to 
[17(1) x U(l) x 1/(1)] x Z 3 CS theory. We expect the 
relation between the Z3 vortices in that theory and the 
Z 3 twist defects in the U(l) x C7(l) x C7(l) theory to be 
similar to the case of the two-component Abelian states. 



VIII. DISCUSSION 

In this paper, we have studied the possibility of ex- 
trinsic twist defects, associated with a symmetry of the 
topological quantum numbers of a topologically ordered 
state. We have explicitly studied several different exam- 
ples, including various kinds of Z 2 and Z 3 twist defects 
in Abelian and non- Abelian topological states. In all of 
the cases we have studied, we found that the projective 
braiding statistics of the twist defects are related to adi- 
abatic modular transformations of a topological state on 
a high genus surface, and we have developed several dif- 
ferent ways of explicitly computing the braiding of the 
twist defects. The close and precise relation of the twist 
defects to topological states on a high genus surface leads 
us to refer to them as genons. 

Most importantly, the braiding statistics of the genons 
are projective, meaning that they are only well-defined 
up to an overall phase. The reason for this is that the 
genons are not intrinsic dynamical quasiparticles of the 
system. This fact allows the possibility of fundamen- 
tally novel behavior, because the projective non-abelian 
braiding statistics are not subject to the same stringent 
constraints as true non-abelian quasiparticles. In this pa- 
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per, we have found several examples of this. The simplest 
case of Z<z twist defects in Abelian states have quantum 
dimension y/n for integer n. For n = 2, these defects are 
Majorana fermions and their braiding is the same as the 
braiding of Ising anyons, up to an overall phase. However, 
as we discussed in Sec. I VIII the braiding of the confined 
defects for n > 2 is inequivalent to the braiding that is 
obtained when the genons are deconfined, because the 
dimension of the braid matrices are different in the two 
cases. In some cases, we even found that after the sym- 
metry is gauged in the theory so that the genons become 
deconfined, it appears that their quantum dimension can 
also change. 

While the non-abelian genons in Abelian states have 
been argued to be insufficient for universal TQC ) 2270 we 
have found that Z2 genons introduced into non-abelian 
states can be used for universal TQC, even when the host 
non-Abelian state is non-universal for TQC. 

In this work, the twist defects are point defects in the 
sense that far away from the defects, no local operation 
can distinguish the presence of the defect, because the 
defects are associated with symmetries of the topologi- 
cal order. It is possible to consider a more general class 
of defects, which are different kinds of domain walls be- 
tween different topological states. Then, junctions be- 
tween different kinds of domain walls may induce some 
non-trivial topological degeneracy and cannot be viewed 
as point defects at all. We leave the detailed study of 
such more general possibilities for future work. 
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Appendix A: Ground state degeneracy with Z-$ twist 
defects via minima of tunneling operators 



Similar to the Zi case, we shall obtain the ground state 
degeneracy of the system by considering the degenerate 
minima of the tunneling operators in (1791) . Since the 
Klein factors are defined differently with to — I even and 
odd, we will separate the discussion into two pieces. We 
assume g < throughout the discussion. 



a. m — I even 

In the absence of twist defects, the whole 1+1D edge 
is effectively in the region A. The operators Ki, K2 and 
/C3, which have eigenvalues ±1, commute with each other. 
We will project the Hilbert space onto the sector where 
K.\, K,2 and /C3 take specific values, since states with 
different values of £1,2,3 are not topologically degener- 
ate. For our convenience, we choose the sector where 
£1,2,3 = 1- Then the conditions of minimizing 6H tU n in 
are 



m<pi + 14*2 + 14>3 = 2irn±, 
l(f>i + m4>2 + I&3 = 27rn^, 
l(f>i + l(f>2 + mdpi — 2nn^, 



(Al) 



where n a ^_ c 6 Z. These equations can be put into a 
diagonal form: 

(to + 2l)<j) + = 2Tr(nf + n% + n^) = 2np + , 

(to - 00(1,-2,1) = 2n(nf - 2n£ + nf) = 2%p( 1 _ 2A) , 

(to - O0(i,o,-i) = 2tt(^ - n£) = 27rp(i i0 ,-i). (A2) 

Different values that <^/'s are pinned to represent differ- 
ent degenerate ground states. Since fix's are subject to 
the periodicity 4>i ~ 4>i + 27r, we should actually use 
e**+, e**d.-a,i) and e^ 1 - -- 1 ) to label the ground states, 
as shown in Sec. IVI Al From (|A2]l . we obtain (jgTj) : 



e^+\p) 



2*ip + 

e #<i,-2,i)|^ = e — —1 — 

6*^(1,0,-1) 1^ = e 

We find that e l< ^+ can take to + 21 different values, while 
each of e^d.-a.i) and e^t 1 - - 1 ) has to— I different choices. 
Thus, in the absence of twist defects, the ground state 
degeneracy of the system on a torus is \(m + 2l)(m— l) 2 \, 
as expected from the U(l) x £7(1) x U(l) CS theory. 

Now we consider the case with n triplets of twist de- 
fects. First of all, by the same reasoning above, we can 
project the Hilbert space onto the sector where ICi, K,\ 
and KJ[ take specific values, since they commute with 
each other. We choose the sector where JCj = JC'j = 
JC'j = 1. The ground states are labeled by the maximally 
commuting set of gauge invariant operators: 

{ e i <t>+( x A i ) e #(i,-2,i)0^i) 6*0(1,0,-1)0*^) e i<t>+(.x Bi ) 
e i<t> + (x Ci ) ^ 6*0(1, o,-i)( a: Bi)+ i 0(-i,i,o)( x c j )l (^3) 

For all operators in region A, the analysis above fol- 
lows straightforwardly leading to the conclusion that each 
of e l ^ +l ^ XA ^ can be pinned to one of to + 21 different 
values, while each of e^d.o.-i)^) an d ^1,-2,1)^) 
can be pinned to one of \m — l\ values. To deter- 
mine the possible values of e ^ +( - XB '\ e^+^oA an( j 
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e #'(i,o,-i)( a: B < )+#' ( '_ 1 , 1 ,o)(^c i ) ) we need t0 write down the 
conditions minimizing SH tun in (|79|) in region B and C: 



27rni 



m<j>'i + Z02 + ^03 
l(pi + m<p2 + ^3 

+ l<j>2 + to0 3 = 2nn 

m<t>" + m + 14>3 

l(f>" + m<t>2 + l<t>3 = 27m 
Z0" + l(j>2 + m<t>3 = 2mi 



B 



<t> fields to different phases which can be written in a 
compact form: 



7T, 



B 



27rnf , 
c 



c 



(A4) 



which can be rewritten in the diagonal form (including 
the position index) 



m + 2l)4>+{x Bi ) = 27r(nf l + rif ' + rif l ), 

m - l)<p' {1 _ 21) (x Bi ) = 2n(nf* - 2n Bi + nf'), 



J.J ( \ i 7 ±f I \ 

TTKpxiXBi) + l<j>2\?Bi) 


+ l<f> 3 (x Bi ) 


= 2-Kn 1 1 - 


1 

1 


-£'i 
2 


7 it ( \ i if I \ 


+ i4> 3 {x Bi ) 


= 27rn 2 - 


1 

1 


-£ 2 
2 


l<j>i{x Bi ) + l<f> 2 (x Bi ) - 


if/ \ 

- m0 3 (x Bi ) 


= 27rnf * - 


1 

V — 


-£ 3 
2 






= 27mf ! H 


1 


-£'/ 




2 


l<f>'{(x Ci ) + m<f>2'(x Ci ) 


+ ^3^0,) 


= 2?™^ H 


1 


-£ 2 ' 

2 




-m^xci) 


= 277713* H 


1 


-£ 3 ' 
9 



(A6) 



m - O0(i,o ,-1)0*^) = 27r(nf i - 
m + 2l)(f> + (xCi) = 27r(ni 4 + n 2 i 
m - O0(i -24) ( x Ci) = 27r(nf s - 
m - 00(1,0,-1) ( x Ci) = 27r(nf s - 



Bi 



), 



The sum of the first three equations in (IA6|) yields 
(m + 2O0+(x Si ) 



2/7 



(A5) 



27r(ni 



3 - fC[ - £ 2 - £ 3 



ir. (A7) 



3 — 1C'^ — JC'2 — /Cg 



Here, n 123 ,n 123 € Z. Thus, we find each of e 
and e i ^+^ xc i l can be pinned to one of to + 21 differ- 
ent values, while e^C 1 ' .- 1 ^ 3 *^"- 1 - 1 . )^ ^ is pinned 
to one of I to — l\ values. This counting gives a total of 
|(m + 2l)(m — Z)| 3 " states. Considering the charge con- 
straint on each defect, which manifests itself in the differ- 
ence between + in different regions (see Sec. I VI A[) . the 

topological ground state degeneracy is |(m+20(m— 3n |- (m - O(0(i -i)( x b,) + 0('_i 1 o){%d)) 



Since, by definition, /C'i/C 2 /C 3 = /C1/C2/C3, 
must be a multiple of 27r in the sector £1,2,3 = 1. Hence, 
we can package the right hand side of this equation above 
into a single integer p^ times 27r. And we can conclude 
that each of e l ^ + ^ XB *' is pinned to one of m + 21 different 
values. The same argument applies for e l *+ ( x °i ) . From 
161), we can also have 



2?r(nf ' - n 



b. m — l odd 

In the absence of the twist defects, the analysis of 
ground state degeneracy is similar to the case above. 
Since the operators £1,2,3 commute with each other, 
we can restrict the system to be in the sector where 
£1.2,3 = 1- Thus, the condition for minimizing 5H tU n 
is the same as ([Ail . In the same fashion, we use the 



4 - JC[ + £ 3 + JC'j - JC'2 1 
2 

(A8) 



Due to the fact that £' 1 £ 2 = — £i£2 and £ 3 £i = 
— JC1JC3, we have K,[ + JC 2 ' — and £ 3 + £" = in the 
sector where £1,2,3 = 1. So we can rewrite the equation 
as 



operators e 



i<t>-\ 



and e 4 "^ 10 '- 1 ' to distinguish the 



(m - (0(i,o,-i) fast ) + ^(-l.i.o)^)) 



degenerate ground states. The ground state degeneracy 
on torus without twist defects is, therefore, given by the 
same formula |(m -I- 21) (m — l) 2 \. 

Now we consider the case with n triplets of twist de- 
fects. The treatment will be slightly more subtle than the 
m — I even case. We still make the projection £1.2,3 = 1 
and identify the maximally commuting set of operators 
as (|A3j) . Thus, all the analysis within region A stays 
the same. Each of e l ^+ ( XA * ' can be pinned to one of 
m + 21 different values, while each of e 1 ^ 10 '- 1 '^- 4 ^ and 
e »<A(i,-2,i)(o^i ( ) i s pinned to one of \m— 1\ values. However, 
in region B and C, £'1 2 3 and £'/ 2 3 cannot be fixed be- 
cause they anti-commute with £1.2,3- With different IC'j 
and £" values, the cosine potentials in ([75)1 will pin the 



27r(nf 



Bi 



,Bi 



1). 



(A9) 



Therefore, e^.a.-o^^-^d.i.o)^) [ s phmed to one 
of I to — l\ different values, just like that in the case with 
to — I even. Then, this counting will give us in total 
of |(m + 20 (to — 0| 3n states. Considering the charge 
constraints on each defect, we conclude the topological 
ground state degeneracy is \(m + 2l)(m — 3 "l- 



Appendix B: Identities of Wilson loop operators 

The definition of the Wilson loop operators are given 
in ([9H1) . Now we will prove (pJTj) . To simplify notation, 
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we consider the case with k = in (|90[) . The Wilson loop 
operators are written as 

W(ai) = ai^273, W{bx) = 73^4^5, 
W(a 2 ) = ^i72«3, W^(6 2 ) = a 3 /3 4 ~f 5 , 
W (a 3 ) = -e^4a 3 , W(b 3 ) = -e&* (Bl) 

Expressed in terms of boson fields, we rewrite W{a\) as 

W(ai) = e #i(^i) e -^i(^2) (B2) 

With the assumption that the charge on each dislocation 
is 0, 0+ is pinned to the one single value globally. From 
(IA1[) . ,we can obtain 



(to - 001 (aMi) - (m - 001 = S^nf 1 - 2 ) 



(B3) 



Thus, 



The proof for W(bi) m ~ l = W(a 2 ) m ~ l = W(b 2 ) m - 1 = 1 
is parallel. Now, we will focus the discussion on W(a 3 ). 
First, we write 

Notice that the 4 operators in this expression commute 
with each other. Thus, 



(a\a 3 ) m ~ l =(— l)e _i ^ T " _ ^^( a: B 1 ) e -i(m-0<;i'i(3;c 1 ) 



X e ~ i ( m ~')<£2( :E A 2 ) e i(TO-()^ + 



(B6) 



Rewriting the second formula in (jAll) as (to— 002+^0+ 
2^71^ , we have 



-i(m — l)4>2(xA 2 ) — e il<t> 



(B7) 



For the operator e ~Km~i)^(x Bl ) e -i{ m -i)4,'{(x Cl ) ^ the 
even/odd-ness of to — Z matters. For m—l even, we obtain 
from $K4$ that (m-O(0 3 (^fli)+0i (a; Cl ))+2Z0 + = 2nZ, 
where Z means an integer, which leads to 

e -K'' n -l)'t>'3( x B 1 ) e -i(m-l)4>'l(xc- l ) _ e i2l<p+ (B8) 

For to — Z odd, 

/C 3 + /C'/ 



(to - O(03(z Bl ) + 0i'(»ci)) + 2Z0+ = 2ttZ + 7T - 



2 

(B9) 



Since JC' 3 JC" = — /C3/C1 = —1, we have K,' 3 + KJ{ = 0. So, 
in this case, 



-i(m-l)4>' 3 (x Bl ) -i(m-l)4>"(x Cl ) _ _ p i2lc/>- i 



(B10) 



Thus, for general m — l, 

e -i(m-l)^' a (x Bl ) e -i(m—l)^"(xci) _ (_ ^ m ~ l e i2l <i>+ _ (Bll) 

Now, we have 

(4a 3 ) m - ! = (_i)»»-'+i e <(m+20* + = (-I)™"'" 1 . 

(B12) 

Considering the phase factor in the definition of W(a 3 ), 
we have proven that W(a 3 ) m ~ l = 1. And, by cyclic 
symmetry of the layer index, we also have W(b 3 ) m ~ l = 1. 
Besides the the Wilson loop operators defined in , we 
will derive similar identities for the Wilson loop operators 
oi 3 Pi^Q and /3374CK6 which is useful in the ID protocol. 
Written in terms of the boson fields, 

«3/3476 = e** 2 '^ 01 ^ - *^^^)^" (^2)^**1(^3) 

— e'-T^T e -iK( x B 2 ) e i(t> 2 ' (xcx ) +i<j>" (xc 2 ) e -#i ( x a 3 ) _ 

(B13) 

Thus, by taking the to — ith power of this expression 
and rearranging the operators while taking all possible 
commutator in account, we have 

(a 3 f3 4l6 ) m - L 

= ( — l) m_ 'e _ '^ m_ '-'*2( :,: B2) e i ( m ~')(^2( a: Ci)+0i'( a; C2)) 
= ( — l) m_ ' e _i ( m ~')*2( :I: B 2 ) e i(m-0(^2'( a; C2)+</ , i'( :I; C 2 )) 

x e - i ( m -O0i( a; A3) e i(ni-O(02'( ;I; C 1 )-02'( a; C 2 )) ; 
= ( — l) m_ 'e _i ( m_ ')'*2( ;1: B 2 ) e -i(m-()^3 , ( a: C2) 

x g-«(« l -O0i( a; A3) e j(» l -O(02'( :I: Ci)-02'( :I: C 2 )) e i(m-O0+ 

(B14) 

In this derivation, we have used the fact that 
e i(m-/)(0 2 (3; Cl )-0 2 Oc 2 )) commute with e - ^" 1- '^^*^) 
e -i(m-l)M*A 3 ) and e *(™-O0+. From (JXTJ or ggj' 
we have e^-W^'^i)-^'^)) = L And fol _ 
lowing a parallel discussion of 1^(03), we obtain 

e -i(m-l)4>' 2 (xB 2 ) e -i(m-l)4>3(xc 2 ) — (_\\m-l £ i2l<j> + anc | 

e _i(m-0MsA 3 ) = e a ^+. Therefore, we conclude 

(a 3 p4-f 6 ) m ~ l = e i{ - m+2l ^+ = 1. (B15) 

By cyclic symmetry of layers, we also obtain 
(p 3l4 a 6 ) m - 1 = 1. 
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Projective representations of the permutation group were 
proposed for arbitrary dimensions in Ref. |57 |. and later 
shown to be inconsistent with locality in Ref. |58| . A differ- 
ent possibility, called projective ribbon permutation statis- 
tics, was found to occur 3+1 dimensions . 11 ' 12 In this pa- 
per, we instead concentrate on projective representations 
of the braid group, which has received little attention in 
the physics literature. 

We would like to clarify that this name has a different 
meaning from the concept of geon in the general relativity 
literature, which refers to a certain type of gravitational 
electromagnetic soliton solution.— However, we also note 
that the topological geon studied in the quantum gravity 
literature^ has some similarity to our genon, since both 
are related to manifold topology and the mapping class 
group. 

For Zi twist defects, we expect this can occur physically 
by proliferating double-defects. 

Strictly speaking, a UMTC describes a topologically or- 
dered state of bosons. The mathematical framework for 
topologically ordered states of fermions is slightly differ- 
ent, as the spins Qi are only defined modulo n. 
To be more precise, the condition above is not in a gauge 
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invariant form and we should require it to be satisfied in a 
certain gauge choice. 

Strictly speaking, <j>Li is left-moving only if the eigenvalues 
of K are all positive. Otherwise it is either right-moving or 
a combination of right- and left-moving modes. Neverthe- 
less, we will continue to use the subscripts L or R even in 
these more general cases. 



See eg Ref. |6l| for a recent derivation. 

This expectation is based on previously known examples 
where Zi gauge theory emerges by proliferating double vor- 
tices in a superfluid. 
Zhenghan Wang, private communication 



